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NOMENCLATURE 


Wing  aspect  ratio 

Wing  chord,  c  =  Wing  mean  chord 


Section  profile  drag  coefficient 

-n-r  A 


Wing  lift  coefficient  - 


for  elliptic  load 


Convected  flow  work  per  unit  area 
Vorticity  moment  of  inertia 


Pressure 


Dynamic  pressure  = 


:u+u)a  +  v8 


Mean  sum  of  the  squares  of  the  turbulent  velocity 
fluctuations 


Radial  distance  from  vortex  axis 


Core  radius 


Wing  half  span 


Time 


U 


Axial  Incremental  velocity  above  free  stream  value 
Free  stream  velocity  (Flight  Speed) 
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v  Rotational  speed 

v  Potential  flow  rotational  speed 

P 

vc  Rotational  speed  at  outer  core  edge 

x  Axial  distance  measured  downstream 

y  Lateral  distance  from  the  plane  of  symmetry 

Greek  Letters 

P  Unknown  constant 

Y  Wing  shed  vorticity 

r  Circulation 

T  Wing  root  circulation 

o 

60^  Fraction  of  induced  drag  coefficient 

i 

AH  Total  head  loss 

Ap  Pressure  increment  above  ambient  pressure 

p  Fluid  density 
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ABSTRACT 

The  effect  of  wing  span  loading  on  the  development  of  fully 
rolled  up  wing  trailing  vortices  is  discussed.  It  is  shown  that 
parabolic  wing  loadings  produce  potential  flow  maximum  core  ro¬ 
tary  speeds  which  are  finite  and  less  than  fifty  percent  greater 
than  the  downwash  speeds  at  the  plane  of  symmetry.  The  develop¬ 
ment  of  turbulent  cores  is  analyzed  and  core  growth  is  predicted 
to  occur  as  the  two  thirds  power  of  time  whereas  the  peak  ve¬ 
locities  fall  off  as  the  inverse  one  third  power.  Axial  flow 
effects  of  the  wing  profile  drag  and  lifting  system  are  shown  to 
lead  to  axial  jets  on  the  vortex  axis  which  may  either  follow  the 
aircraft  or  exceed  the  free  stream  velocity  depending  on  the 
ratio  of  profile  drag  to  induced  drag. 
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INTRODL'CTION 

The  twin  tornado-like  vortex  systems  created  behind  a 
lifting  airplane  wing  have  recently  attracted  much  attention 
because  of  the  potential  hazard  to  aircraft  that  encounter 
them.  The  primary  danger  appears  to  be  the  development  of 
high  rolling  rates  **or  following  aircraft  that  may  fly  di¬ 
rectly  into  the  center  of  the  high  velocity  rotating  air  mass. 
Transverse  flight  through  the  whirling  vortex  system  may  also 
cause  high  gust  loads  on  the  encountering  aircraft.  Of  par¬ 
ticular  concern  are  the  upsetting  effects  on  aircraft  landing 
at  busy  air  terminals  where  sufficient  maneuvering  altitude 
may  not  be  available  to  recover  from  the  rolling  or  pitching 
motions  induced  by  the  vortices.  Presently  several  agencies 
of  the  government  are  studying  the  problem  seeking  understand¬ 
ing  of  the  basic  aerodynamics,  means  for  detecting  the  presence 
of  the  vortices  and  also  means  for  reducing  their  intensity. 

There  is  currently  some  debate  on  the  longevity  of  the  vortices 
once  produced  and  various  studies  have  been  undertaken  to 
measure  and  predict  the  vortex  decay  rates.  It  appears  that 
atmospheric  conditions,  in  particular  the  turbulence  level, 
have  much  to  do  with  the  longevity. 

The  present  paper  is  concerned  with  the  basic  aerodynamic 
patterns  which  govern  the  development  and  decay  of  the  vortex 
system  downstream  of  the  generating  aircraft.  Calculations  are 
carried  out  of  the  vortex  downwash  field  for  the  case  of  para¬ 
bolic  wing  loading  and  comparison  is  mrde  with  that  of  e.lliptically 
loaded  wings.  For  elliptic  load  conditions  the  effect  of  the 
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wing  profile  drag  on  the  axial  flows  is  considered  and  finally 
the  conditions  are  derived  under  which  self  similar  wake 
development  can  occur. 

GENERAL  DISCUSSION 

Prandtl,  Kaden  and  Betz,  References  1,  2  and  3> 
respectively  described  the  potential  flow  to  be  expected 
behind  a  lifting  three  dimensional  wing.  Betz,  however,  pre¬ 
sented  the  first  calculation  method  which  described  quantita¬ 
tively  the  velocity  and  vorticity  fields  far  behind  an 
elliptically  span  loaded  wing.  The  vortex  sheet  shed  from 
the  wing  trailing  edge  was  ultimately  shown  to  roll  up  into  a 
pair  of  spiral  sheets  with  an  infinite  number  of  convolutions 
so  that  the  vorticity  and  tangential  velocity  were  essentially 
continuous  functions  of  the  distance  from  the  center  of  the 
Spiral.  The  vorticity  was  computed  to  be  well  spread  out  so 

that  the  singular  tangential  velocity  at  the  vortex  system 

-A  _ 

Center  was  found  to  vary  as  r  2  rather  than  r  as  would  be 
the  case  for  a  concentrated  line  vortex. 

The  modification  of  the  potential  flow  by  viscosity  or 
turbulent  mixing  was  not  discussed  till  1954  when  Squire 
Reference  4  attempted  to  calculate  the  decay  rate  of  the 
vortex  system  as  a  function  of  time  after  generation.  Squire’s 
assumption  of  an  initial  concentrated  line  vortex  and  other 
simplifying  assumptions  such  as  the  constancy  of  the  turbulent 
eddy  viscosity  make  his  model  somewhat  questionable.  Batchelor, 
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Reference  5,  has  presented  an  excellent  description  of  the 
rolling  up  process  including  effects  of  the  wing  viscous  wake 
however-,  his  very  enlightening  work  on  the  development  of  the 
central  core  of  the  vortex  assumes  laminar  flow  and  again  he 
starts  with  the  velocity  field  of  a  concentrated  line  vortex. 

The  vortex  wakes  are  generally  found  to  grow  at  rates  much 
slower  than  those  of  non-rotating  wakes  and  it  is  now  known 
that  the  rotation  produces  a  certain  stabilizing  effect  on  the 
system  and  appears  to  cause  a  reduction  in  turbulent  eddy 
diffusivity  at  least  in  the  radial  direction.  Most  authors 
dealing  with  the  vortex  wake  problem  have  taken  a  concentrated 
line  vortex  as  a  starting  point.  Clearly  the  initial  conditions 
of  energy  distribution  and  field  stability  must  play  important 
roles  In  the  subsequent  system  decay  and  it  is  therefore  im¬ 
portant  to  establish  these  conditions  and  to  investigate  their 
consequences.  We  begin  by  looking  into- the  effect  of  altered 
span  loading  on  the  final  potential  flow  vortex  velocity  field. 

Effect  of  Span  Loading 

The  procedure  of  Betz  has  been  applied  to  the  case  of 
elliptic  loading  by  Betz  himself.  Reference  3  and  to  the  case 
of  a  rectangular  wing  by  Mason,  Reference  6.  In  this  method 
it  is  assumed  that  the  development  of  the  vortex  spiral  pro¬ 
ceeds  as  a  two-dimensional  flow  that  is,  the  fore  and  aft 
influence  at  a  downstream  station  is  shown  to  be  small  and  is 
neglected.  The  lateral  position  of  the  center  of  gravity  of 
the  vorticity  for  each  side  of  the  vortex  system  i6  shown  to 
be  constant  as  is  the  moment  of  inertia  of  vorticity  about 
t.ne  lateral  center  of  gravity.  These  considerations  are  shown 


V  s«i 


fi  I^^2t&£saV?**fg*£ 


HYDRONAUTICS,  Incorporated 

-5- 

to  follow  from  Kelvin's  original  theorems.  At  this  point  Betz 
introduces  the  approximately  correct  assumption  that  the  moment 
of  inertia  of  each  vorticity  group  about  its  own  center  of 
gravity  is  constant  during  the  spiral  roll  up.  Specifically, 
the  vorticity  shed  from  a  region  of  the  wing  span  from  the 
tip  to  any  inboard  station,  y,  rolls  up  into  a  spiral  having 
a  radial  moment  of  inertia  of  vorticity  about  the  group  center 
of  gravity  equal  to  the  moment  of  inertia  of  the  original 
vorticity  distribution  shed  from  the  wing.  Thus  if  the  in¬ 
itial  span  loading  is  known  the  potential  flow  in  the  rolled 
up  spiral  downstream  at  infinity  can  be  computed. 

In  rjost  cases  the  shed  vortex  distribution  is  singular 
at  the  wing  tip  having  an  (s-y)  2  variation  at  that  point, 
here  s  is  the  half  span  distance  and  y  is  the  distance  to  a 
point  on  the  span  measured  from  the  center  line  or  wing  root. 
The  calculations  of  the  rolled  up  spiral  velocity  field  show 
corresponding  singular  behavior  in  the  core  or  center  of  the 
spiral.  This  effect  will  be  discussed  in  more  detail  in  the 
following  sections.  It  is  instructive  to  calculate  the  rolled 
up  field  for  the  case  of  parabolic  loading  where  there  is  no 
singular  behavior  at  the  wing  tip.  We  therefore  take  the  span 
loading  as 

r(y)  *  rQ  (1  -  §■)  (i) 

where  r  is  the  local  circulation,  r»  is  the  wing  root  cir¬ 
culation,  y  is  the  lateral  distance  from  the  root  and  s  is  the 
semi  span.  Hie  vorticity,  y ,  shed  in  the  wake  close  behind  th.- 
wing  is  now 
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Pr  y 
o 


(2) 


Parting  the  subgroup  of  vorticity  shed  from  y  to  s  along  the 

spur,  we  compute  the  location  of  the  center  of  gravity,  y 

eg 

■Arid  the  moment  of  inertia,  J,  as 


5 

/ 


Cg  =  r[yT  /  'Y(y,)y,dy,  = 


2s  ‘ 

I-* 

i 

> 

W 

- - > 

3 

l  -  (y/s)3] 

and 


J(y) 


eg 


c 

■/ 


sar 

Y(y'-ycg)ady'  =  — 


(*) 


Tr;  tne  rolled  up  spiral,  this  group  of  Cortices  must  be  contained 
within  a  radius,  r,  and  have  the  same  moment  of  inertia  about  the 
core  center,  therefore 


Vi.i )  =  r(r)  and  J(y)  =  J(r) 


!,,rom  (h)  and  (l)  we  may  write 


dr.  dr, 
dF*  ‘  to?*  ■ 


SV  . 

-p-ay 


(5) 


(6) 


arid 
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dj j  dJ (r)  , 
3— dy  =  -r*-  ‘  dr 
dy  17  dr 


(7) 


By  definition 


J(r) 


r 

■/ 


rsdr 

dr 


(8) 


and  by  combining  (4),  (6),  (7)  and  (8)  we  obtain  a  relation¬ 
ship  between  r  and  y  as 


Pa  ya 

7*'  Pr  + 


»[■ 


(1  +  *>• 


(9) 


This  functional  relationship  is  plotted  in  Figure  1  together 
with  the  corresponding  functions  for  the  elliptically  loaded 
wing.  Note  that  r/s  approaches  2/3  for  y/s  =  Oj  this  value 
corresponds  to  the  center  line  of  the  double  vortex  system 
for  the  case  of  parabolic  span  loading.  As  pointed  out  by 
Bet/,  the  interaction  of  the  two  systems  which  has  been  ne¬ 
glected  would  undoubtedly  somewhat  modify  these  results  for 
large  y/s. 

The  circulation  distribution  in  the  rolled  up  spiral 
vortex  can  be  obtained  as  a  function  of  r/s  as 

r(r)  =  ro  U  - 


;4,a  daafaaliaktStltMii  at  maCi3MKI 


WP?!W^W. 
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Mear  the  axis  (9)  can  be  approximated  as 


<5>  ■  a  -  H) 


(ii) 


hence  the  circulation  near  the  core  is  approximately 

r 


z  r .  w 

**(8)  -  s 


(12) 


The  complete  function  is  plotted  in  Figure  1. 

As  a  result  of  the  finite  slope  of  the  circulation 
distribution  at  r  =  0  we  infer  that  the  velocities  at  the 
vortex  center  line  would  be  finite.  The  rotary  velocity 
field  is  given  simply  for  axisymetrical  flows  as 


v  =  ri£l 

2m 


(13) 


and  at  the  axis  the  rotary  speed  is  thus 

2l\ 


ruO 


_ o 

7TS 


(14) 


The  complete  downwash  pattern  in  the  plane  of  the  vortex  pair 
has  been  calculated  approximately  from  the  above  relations  and 
is  plotted  in  Figure  2;  also  shown  is  the  downwash  computed  for 
the  elliptic  span  load  distribution.  It  should  be  noted  that 
the  center  of  gravity  of  the  vortex  system,  (taken  as  the  center 
of  the  vortex  spiral)  for  the  parabolic  span  loading  is  located 
at  the  two-thirds  span  position  whereas  that  for  the  elliptical 
loading  is  at  a  point  tt/4  times  the  semi  span. 
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Pigure  2  brings  out  a  major  conclusion  of  this  work, 
namely  the  maximum  core  velocity  for  parabolic  wing  loading  is 
finite  and  less  than  50#  higher  than  the  downwash  velocity  at 
the  centerline.  If  one  considers  an  aircraft  encountering  these 
velocity  distributions,  that  is,  flying  with  its  fuselage  along 
the  axis  of  the  vortex  spiral  it  appears  apparent,  nonetheless, 
that  the  induced  rolling  moment  would  be  greater  for  parabolic 
loading  than  for  elliptic.  The  elimination  of  the  high  rotary 
speeds  on  the  axis  is  consequently  of  little  help  in  alleviating 
the  flight  hazard.  It  becomes  clear  that  the  reduction  in 
rolling  moments  on  encountering  aircraft  must  come  from  reduc¬ 
tions  in  the  rotary  speeds  well  out  from  the  vortex  axis.  Far 
downstream  from  the  aircraft  the  vortex  system  can  be  expected 
to  decay  under  the  influence  of  viscosity  or  turbulent  diffu¬ 
sion  of  vorticity  and  it  is  important  to  determine  whether  the 
velocity  field  is  substantially  altered  only  near  the  axis  in 
a  central  core  or  whether  the  angular  momentum  of  more  distant 
regions  can  be  reduced.  The  answer  to  this  question  involves 
the  stability  of  the  system  and  some  consideration  of  these 
questions  will  be  given  in  the  following  sections. 

Returning  to  the  elliptical  load  distribution  case,  it  is 
seen  in  Figure  2  that  very  high  rotational  speeds  are  predicted 
close  to  the  axis.  Confirmation  of  the  velocity  field  calcula¬ 
tions  has  been  obtained  in  several  places.  References  6  and  7» 
and  recently  tests  sponsored  by  the  Langley  Research  Center  of 
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the  NASA  have  provided  the  data  shown  in  Figure  3*  These  hot 
film  probe  data  were  obtained  for  a  model  of  the  Boeing  747  air 
craft  at  its  cruise  flight  condition.  The  tests  were  carried 
out  in  the  HYDRONAUTICS,  Incorporated  Ship  Model  Basin  at  a 
Reynolds  number  of  1.5  x  10®  based  on  wing  mean  chord.  Note 
that  at  a  distance  of  63  wing  spans  behind  the  aircraft  cor¬ 
responding  to  2.3  miles  full  scale,  a  peak  velocity  of  over 
six  times  the  centerline  downwash  value  was  measured  at  a 
radius  of  just  over  one  percent  of  the  wing  span.  Although 
there  is  scatter  in  the  data  they  yield  a  remarkably  good  fit 
to  the  theoretically  derived  curve.  No  substantial  reductions 
in  outerfield  rotary  speed  is  noted  at  the  110  span  (4.1  miles) 
downstream  position. 
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Conslderations  of  Axial  Flows 

The  type  of  radial  circulation  distribution  associated  with 
potential  flow  in  a  rolled  up  spiral  sheet  is  generally  a  mono¬ 
tonic  increasing  function  of  radial  distance  from  the  spiral 
center  (see  Figure  1).  Such  distributions  in  the  absence  of 
axial  flows  are  known  to  be  stable  with  respect  to  radial  dis¬ 
turbances  (Rayleigh,  Reference  8  and  Synge,  Reference  9) •  One 
may  ask  how  then  a  turbulent  core  develops  and  grows  with  down¬ 
stream  distance.  One  answer  can  be  found  in  the  presence  of  the 
axial  flows  associated  with  the  profile  drag  of  the  wing  in 
combination  with  the  axial  flows  produced  by  the  rotary  flow 
itself.  The  latter  flow  has  been  discussed  by  Batchelor, 
Reference  5*  who  derived  the  equations  for  the  axial  flow 
associated  with  rotary  flow  fields. 

The  pressure  equations  of  pertinence  of  axisymmetric 
rotary  flows  are: 


Ap  =  -pUu 


AH 


(20) 


and 


Ap  -  - 


dr  +  constant 


(21) 


r 

in  which  Ap  is  the  pressure  increment  above  ambient  pressure,  p 
is  the  fluid  density  and  AH  is  the  total  head  loss.  The  constant 
in  (21)  is  zero  to  first  order.  If  we  assume  the  head  loss,  AH, 
occurs  as  a  result  of  the  profile  drag  losses  on  the  wing  it  is 
possible  to  relate  AH  and  profile  drag  through  the  axial  momentum 
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equation  and  the  assumption  that  the  wing  wake  rolls  up  with 
the  vortex  system  so  that  the  axial  force  distributes  itself 
exactly  as  does  the  vorticity.  This  assumption  is  a  rather 
natural  one  because  the  trailing  vortex  lines  are  essentially 
embedded  in  the  wing  viscous  wake  as  both  leave  the  wing  trail¬ 
ing  edge.  It  therefore  follows  that  the  low  energy  wake  air 
will  f.-llow  the  axial  vortex  lines  during  the  roll  up  process. 
We  therefore  take  the  profile  drag  on  an  outboard  portion  of 
the  wing  and  assume  that  this  force  appears  in  the  rolled  up 
spiral  within  a  radius  defined  by  the  Betz  procedure,  that  is, 
the  circle  within  which  the  vorticity  shed  from  the  correspond¬ 
ing  portion  of  the  wing  would  appear.  Following  this  argument 
we  write 


+  6C  )qcdy  = 


r 

-2„J 


rAP  +  pUu]rdr 


(22) 


Here  CD  Is  the  wing  local  profile  drag  coefficient,  q  is  the 
P 

stream  dynamic  pressure,  c  is  the  local  wing  chord,  and  6C  is 

i 

some  small  fraction  of  the  lost  momentum  associated  with  lift 
which  appears  in  the  rolled  up  spiral  within  the  circle  of 
radius  r.  Note  that  the  momentum  defect  in  the  spiral  for 
small  r  as  a  result  of  the  drag  due  to  lift  is  only  a  fraction 
of  the  total  induced  drag  because  the  bulk  of  this  drag  is 
distributed  as  kinetic  energy  well  outside  the  cores.  6C 

Di 

should  be  independent  of  profile  drag  and  can  be  evaluated  by 

setting  -  0  and  CD  =0. 

P 


^  -V Jta ;-^/,'cr%.  ,  _  ,.  „ /^^vtOT^prrr-.v. *•  fr 'igm? m3F^*$*GGr**&&^?W& ’W«'»3r 


HYPi-ONAHTICS,  Incorporated 


-13- 


If  we  now  differentiate  (22)  and  make  use  of  Betz's 
approximate  relation  connecting  r  and  y  for  the  case  of 
elliptic  span  loading 


r/s  =  ~  (1  -  y/s) 


(23) 


we  obtain  making  use  of  (20) 


(CDp  ♦  6CDi„c  ■¥[*£  +  *« 


(2M 


From  the  independence  of  profile  drag  and  drag  due  to  lift 

3CC  qc 

*H(r)  .  far  '  I 

Note  that  the  equation  indicates  a  singular  behavior  of  the 
total  head  on  the  axis.  In  reality  of  course,  laminar  or 
turbulent  mixing  would  preclude  such  a  possibility,  however, 
the  rolling  up  process  is  clearly  shown  to  result  in  a  con¬ 
centration  of  low  energy  air  along  the  vortex  axis.  Using 
Betz's  result  for  elliptic  loading,  the  potential  flow 
velocity  field  for  small  r  is  given  by 


o  /2rsx' 
2ir(— ) 


(26) 


hence  (21)  becomes 


AP  =  - 


3ro" 

feirasr 


(27) 


at  ~vi  ^T*jT^v'f  iri 
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Making  use  of  this  with  (25)  and  (20)  enables  the  calculation  of 
the  axial  velocity  distribution  as 


u 

U 


3r  2 
o 


l6vau2s2 


3Cd  c 


8  vs 


s_ 

r 


(28) 


Expressed  in  terms  of  the  usual  parameters ,  the  axial  velocity 
distribution  in  the  rolled  up  core  is  found  to  be 


u 

U 


jV 

|_4v*a2 


3Cd  c ( r/s ) 3 

_ P 

8irAc 


r 


(29) 


where  A  is  the  wing  aspect  ratio,  C,  the  wing  lift  coefficient, 
and  c/c  the  ratio  of  local  wing  chord  to  mean  chord,  c. 

It  can  be  seen  from  Equation  (29)  that  the  incremental 
axial  velocity  can  be  either  forward  or  rearward  depending  on 
the  lift  and  profile  drag  coefficients  of  the  wing.  Also  ap¬ 
parent  from  (29)  is  the  tendency  for  the  axial  flow  to  focus 
near  the  vortex  axis  as  indicated  by  the  r-1  dependence. 

It  should  be  noted  that  Equation  29  should  be  expected  to 
hold  only  at  relatively  short  distances  behind  the  aircraft  where 
the  spiral  system  is  well  rolled-up  but  before  the  field  has  been 
greatly  spread  and  reduced  in  intensity  by  turbulent  diffusion. 
Much  farther  back  the  processes  described  by  Batchelor,  Reference 
5,  must  come  into  play  and  the  field  would  always  present  a  wake 
following  the  aircraft.  That  is,  both  the  profile  drag  and  drag 
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due  to  lift  must  ultimately  show  up  downstream  as  a  momentum 
defect,  hence,  a  following  wake  must  occur  even  if  the  system  is 
broken  up  by  three  dimensional  instabilities  of  the  vortex  pair. 


Stability  of  the  Rotating  System 


Having  obtained  an  approximate  description  of  the  flow  fields 
to  be  expected  in  the  rolled-up  trailing  vortex  system  at  distances 
of  a  few  span  lengths  behind  the  aircraft,  it  is  now  possible  to 
discuss  the  stability  of  the  flow.  It  is  of  interest  to  determine 
whether  the  flow  field  or  certain  portions  of  the  flow  field  have 
a  natural  tendency  to  produce  turbulence  by  feeding  energy  from 
the  mean  flow  motion  into  random  motion  or  whether  the  system  is 
one  which  tends  to  damp  out  turbulent  fluctuations.  Unfortunately., 
there  is  some  disagreement  among  investigators  as  to  general  con¬ 
ditions  for  stability  of  rotating  flows.  Most  pertinent  to  the 
present  case  is  the  work  of  Ludwieg,  Reference  10,  who  has  pre¬ 
sented  a  general  stability  map  for  rotating  flow  which  involves 
two  parameters,  namely,  the  dimensionless  velocity  gradients  de¬ 
fined  as 


and 


(30) 


(31) 


Ludwieg 's  condition  for  stability  is  given  as 


(1  -  C  )(1  -  C  s) 

x  C 0  CO 


(32) 
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This  relation  is  plotted  in  Figure  4  with  the  stable  and  unstable 
regions  noted.  Also  shown  is  the  stability  boundary  found  in  an 
experimental  study  making  use  of  a  test  arrangement  having  rota¬ 
ting  cylindrical  walls  with  provision  for  axial  translation  of  the 
inner  one.  Introducing  now  the  results  derived  for  the  rolled-up 
vortex  sheet  of  an  elliptically  loaded  wing,  equations  26  and  29> 
there  is  obtained 


C  =  -1/2,  C 
cp  5  z 


3Cls 

3Cd 

p  1 

U] 

2s\ 

i  4tt4  A8 

8irAc  1 

i  2CJ 

3r) 

(33) 


Considering  Figure  4  and  Equation  33  it  can  be  seen  that  the 

flow  field  is  predicted  to  be  unstable  on  the  vortex  axis  (C  =  -<») 

z 

but  to  be  stable  in  the  outerfield  as  C  falls  below  a  value  of 

z 

0.72.  Substitution  of  any  reasonable  values  in  equation  33  shows 
that  the  region  of  instability  is  confined  to  an  extremely  small 
radius,  but  nevertheless  a  region  does  exist  capable  of  initiating 
turbulent  motions  which  could  ultimately  produce  a  turbulent  co^e 
and  modify  the  potential  flow  field. 


It  should  be  noted  that  Howard  and  Gupta,  Reference  11,  have 
criticized  Ludwieg's  result  on  the  grounds  that  the  analysis  is 
not.  mathematically  rigorous.  In  their  own  work  they  failed  to 
find  any  general  condition  for  stability  to  non-axi symmetric  dis¬ 
turbances,  however,  as  they  pointed  out,  this  does  not  imply 
instability  either.  DonaJdson,  Reference  12,  derives  a  stability 
criterion  which  predicts  instability  for  all  values  of  less 
than  -0.37  when  no  axial  flow  is  present  but  his  calculation  also 
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shows  markedly  reduced  Reynold's  stress  for  the  case  under  con¬ 
sideration  of  C  =  -0.5-  Donaldson's  method  is  not  a  rigorous 
cp 

proof  as  it  contains  certain  arbitrary  assumptions  concerning  the 
invariance  of  triple  correlation  coefficients  of  the  fluctuating 
velocity  components  under  changing  conditions  of  the  mean  flow. 

It  appears,  therefore,  that  the  question  of  whether  the  basic  flow 
field  having  a  v  ~  r  2  dependency  is  stable  or  not  must  await 
further  analysis  or  controlled  experiments.  Ludwieg's  experi¬ 
ments  summarized  in  Figure  4  predict  stability  for  the  case  in 
question  and  are  quite  persuasive.  We  will  proceed,  therefore, 
with  some  reservations  on  the  assumption  that  Ludwieg's  theory 
is  correct. 

The  complete  stability  of  the  flow  field  would  mean  that 
all  initial  turbulence  would  decay  and  no  additional  turbulent 
energy  would  be  extracted  from  the  field;  consequently  only 
molecular  diffusion  of  vorticity  would  occur.  However,  because 
of  turbulence  produced  in  the  low  stability  region  of  the  axis 
or  more  probably  because  of  the  initial  turbulent  energy  from 
the  wing  profile  drag  which  is  concentrated  on  the  axis,  a  tur¬ 
bulent  core  develops  which  modifies  the  potential  flow  field  and 
in  some  way  establishes  conditions  of  instability  necessary  for 
the  conversion  ,f  rotary  motion  into  turbulence.  In  the  turbu¬ 
lent  core  region  eddy  viscosity  must  certainly  be  much  larger  than 
molecular  viscosity  even  though  substantial  suppression  of  radial 
eddy  mixing  occurs  as  a  result  of  the  stable  rotary  velocity  field. 
In  the  region  outside  the  core  only  molecular  viscosity  is  avail¬ 
able  to  transmit  shear  stress  and  as  a  result  the  outer  field  must 
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remain  substantially  unaltered  until  it  is  engulfed  by  the  expand¬ 
ing  turbulent  core.  This  situation  '.n  unlike  the  case  of  laminar 
decay  (Reference  13)  in  which  the  ou<er  field  decelerates  the  inner 
field  by  viscous  shearing  stresses  and  in  which  the  internal  angular 
momentum  is  transported  outward  to  infinity  without  local  increases 
in  rotary  speed.  Rather  in  the  present  case  the  original  angular 
momentum  within  the  core  radius  must  be  largely  conserved  so 
thit  some  local  increase  in  rotary  speed  must  be  anticipated 
as  the  turbulent  shearing  stresses  modify  the  potential  flow 
rotary  field.  If  the  turbulent  shearing  stresses  were  large 
or  if  sufficient  time  is  allowed ,  the  core  region  would  tend 
toward  the  minimum  kinetic  energy  state  consistent  with  the 
conservation  of  angular  momentum.  This  condition  is  the  stress 
free  state  of  uniform  rotation  and  is  of  course  only  approach¬ 
able  as  a  limiting  condition  however  it  is  instructive  to 
bring  out  the  consequences  in  the  flow  of  these  tendencies. 

The  velocity  and  circulation  distributions  for  this  limiting 
case  are  shown  as  the  solid  lines  in  Figure  5*  It  should  be 
noted  that  conservation  of  momentum  requires  core  edge  veloci¬ 
ties  to  increase  relative  to  the  original  potential  flow  values 
and  therefore  a  circulation  "overshoot?'  occurs  as  found  in 
References  7  and  8  for  the  case  of  a  line  vortex.  As  the  core 
grows,  it  engulfs  a  portion  of  the  outer  field  having  a  lower 
angular  speed  than  that  of  the  core  hence  the  angular  core 
speed  continually  diminishes  with  time.  The  tendency  toward 
the  circulation  discontinuity  provides  the  necessary  instability 
(dr2/dr  <  0)  to  generate  turbulent  motions  at  the  advancing  core 
boundary,  however  the  calculation  of  the  final  velocity  field  is 
beyond  the  scope  of  the  present  report;  it  is  possible  however 
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that  the  methods  of  References  7  or  8  can  be  applied.  The 
dashed  lines  shown  in  Figure  5  represent  the  author* s  estimate 
of  tne  form  of  the  rotary  velocity  and  circxilation  distribu¬ 
tions.  As  pointed  out  by  Govindaraju  and  Saffman  (Reference  8) 
the  experimental  verification  of  the  velocity  field  will  be 
extremely  difficult  because  of  the  usual  unsteadiness  of  vortex 
flows  which  would  tend  to  smear  out  any  steep  velocity  gradients, 
Finally  it  must  be  remarked  that  the  speed  of  advance  of  the 
turbulent  front  into  the  unperturbed  fluid  must  be  affected  by 
the  rotation  induced  stability  of  that  region  and  one  should 
expect  the  growth  to  be  substantially  less  rapid  than  that  of 
simple  wakes  such  as  those  produced  by  bodies  of  revolution, 
disks,  etc. 


Vortex  Core  Growth 


As  indicated  by  Squire,  Reference  4,  it  is  probable  that 
the  core  grows  in  a  manner  satisfying  certain  self  similar 
flow  relationships  and  it  is  of  interest  to  examine  the  physics 
of  the  flow  for  their  existence.  As  conditions  for  the  self 
similar  development  we  write 


vc  f(r/rc) 


(34) 


AH  = 


PV  g(CD^/CLa,  r/rj 


u  = 


h(CD  /V’  r/rc> 
P 


Ap  = 


pVo‘  t(CDp/CL’’  r/V 


(35) 

( 36) 
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Equations  (36)  and  (37)  follow  from  (31*),  (35)*  (20)  and  (21). 
The  rate  of  turbulent  diffusion  is  known  to  be  proportional 
to  the  root  mean  square  of  the  turbulent  velocity  fluctuations 
and  we  therefore  write 


dr 

c 

~dt 


=  /3 


(37) 


where  /3  is  a  constant  which  should  depend  in  an  inverse  way  on 

the  stability  of  the  outer  velocity  field  and  q3  is  the  mean 

intensity  of  turbulence  in  the  core.  The  turbulent  energy 

arises  from  the  change  in  rotary  kinetic  energy  from  the  initial 

potential  flow  state  co  the  final  state  but  may  be  reduced  by 

dissipation  or  altered  by  energy  flow  convected  by  the  axial 

currents.  Ignoring  both  latter  losses  for  the  moment,  we 

calculate  the  turbulent  intensity  within  the  circle  of  radius 

r  as 
c 


qa 


-  v8  )rdr  ~ 


sr 

c 


(JS) 


O 


in  which  v^  is  the  initial  potential  flow  rotary  speed  and  we 

have  made  use  of  Equations  (26)  and  (3^).  It  is  seen  that  the 

turbulent  Intensity  falls  off  as  the  inverse  of  the  core  radius 

-2 

whereas  for  a  line  vortex  the  intensity  would  vary  as  rc 
Substitution  of  (38)  i  (37)  yields 


dr 

c 

■  '■v# 

dt 


(39) 
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and  upon  Integration  we  obtain 


r 

c 


2/3  S-V3  t2/3 
1  o 


(40) 


We  conclude  from  -,'0)  that  the  conditions  for  self  similar 

P/-3 

growth  are  satisfied  if  the  core  radius  grows  like  t  '  or, 

2/3 

because  t  =»  x/U,  like  x  '  .  The  maximum  rotary  speed  In  the 
core  is  proportional  to  v  and  we  may  therefore  write 

C  r  r  2/3 

'max  ~  Tf==-  ~  “i/3°tl/3  C1* 

V  src  st' 

-1/3 

and  hence  the  maximum  velocity  also  varies  as  x  .  These 

results  are  in  contrast  to  the  findings  for  line  vortices, 

1/2  -1/2 

r  -  x  and  v  ~  x  ,  solely  because  of  the  initial 
c  max 

vorticity  distribution  differences  between  the  two  cases. 

It  is  important  to  note  that  the  result- given  by  Equation  (40) 

should  be  generally  valid  for  real  flows  because  most  vortex 

motions  are  generated  by  rolling  up  of  spiral  vortex  sheets 

-3/2 

in  which  the  initial  vorticity  should  vary  like  r  as  in 

the  present  calculation.  It  should  also  be  mentioned  that 

many  correlations  of  flight  test  data  have  been  attempted 
1/2 

using  the  x  relationship  which  now  appears  to  be  in  error. 
Other  calculations  such  as  the  speed  of  vortex  rings  having 
turbulent  cores  should  also  be  reexamined  in  light  of  the 
above  results. 
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With  dissipation  and  axial  flow  present  the  similarity 
is  not  destroyed  because  (1)  the  dissipation  proceeds  at  a 
rate  proportional  the  basic  energy  release  and  (2)  the  axial 
convection  of  energy  is  constant  to  first  order.  The  dissipa¬ 
tion  in  a  turbulent  flow  can  be  generally  written  (Reference  13) 


. 

at 


t  ~ a \ 3/2 


(42) 


and  upon  integration  introducing  Equation  (40)  we  compute  the 
energy  dissipated  by  viscosity 


Aq8/2  ~  - 


(*3) 


Comparison  of  (43)  with  (38)  shows  that  the  dissipation 
is  directly  proportional  to  the  kinetic  energy  released  and 
therefore  the  self  similar  development  is  preserved. 

The  convection  of  turbulence  and  flow  work  by  the  axial 


currents  can  be  written 


9tl  ax  )  r 


{♦  /'■ 


(q3  )rdr 


(44) 
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lntroducing  Equations  (36),  (37)  and  (38)  into  the  above  we 
find  that  the  terms  within  the  braces  become  independent  of 
r  or  in  physical  terms  the  lateral  convection  of  energy  in  the 
core  is  constant  to  first  order.  As  a  consequence  the  lateral 
energy  flow  does  not  affect  the  energy  balance  and  the  conditions 
for  self-similar  development  are  not  violated. 

CONCLUDING  REMARKS 

The  development  of  turbulent  vortex  cores  In  a  twin  spiral 
rolled  up  vortex  sheet  behind  a  lifting  wing  has  been  discussed. 
The  effect  of  wing  span  loading  is  shown  to  be  important  in  re¬ 
gard  to  the  peak  rotary  speeds  developed  in  the  vortex  Bystem. 
Calculations  indicate  that  parabolic  span  loading  on  the  wing 
produces  core  maximum  rotary  speeds  less  than  fifty  percent 
higher  than  the  centerline  downwash  velocity. 

Under  the  plausible  assumption  that  the  viscous  wing  wake 
rolls  up  in  the  vortex  spiral  with  the  trailing  vortices,  it 
is  found  that  in  the  early  stage  of  vortex  development  axial 
currents  are  formed  at  the  vortex  center  which  may  be  either 
greater  or  less  than  the  flight  speed  depending  on  the  ratio  of 
profile  drag  to  induced  drag. 

For  aircraft  having  elliptic  span  loadings  or  indeed  for 
all  vortical  flows  where  the  shed  vorticity  has  an  inverse 
square  root  type  singularity  at  the  generating  surface,  the 
turbulent  core  of  the  rolled-up  vortex  system  is  pxedicted  to 
grow  as  the  two-thirds  power  of  time  or  the  downstream  distance. 
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The  peak  rotary  speeds  fall  off  as  the  inverse  one  third  power. 
These  results  are  believed  to  be  new  and  in  contradiction  to 
earlier  analyses  made  assuming  the  vortex  system  originated  as 
discrete  concentrated  vortices. 

Finally,  it  is  concluded  that  the  rolling  moment  induced 
on  aircraft  entering  the  vortex  wake  should  be  little  affected 
by  the  core  decay  until  the  core  diameter  reaches  a  size  com¬ 
parable  with  the  encountering  aircraft's  wing  span. 
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RADIAL  DISTANCE  OUT  FROM  THE  CORE  CENTER  -  HALF  SPANS 


FIGURE  1  -  CIRCULATION  DISTRIBUTION  FOR  ELLIPTIC  AND  PARABOLIC 
LOADINGS 
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FIGURE  3  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  WAKF.  DOWNWASH 


FIGURE  4  -  LUDWIEG'S  STABILITY  MAP 
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figure  5  -  LIMITING  VELOCITY  AND  CIRCULATION  PROFILES 


